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GRB 021211 as a Faint Analogue of GRB 990123: Exploring the 
Similarities and Differences in the Optical Afterglows 1 

Stephen T. Holland 2,3,4 , David Bersier 5 , J. S. Bloom 6,7 , Peter M. Garnavich 4 , Nelson 
Caldwell 6 , Peter Challis 6 , Robert Kirshner 6 , Kevin Luhman 6 , Brian McLeod 6 , & 

K. Z. Stanek 6 

ABSTRACT 

We present BVRcJHK s photometry of the optical afterglow of the gamma- 
ray burst GRB 021211 taken at the Magellan, MMT, and WTYN observatories 
between 0.7 and 50 days after the burst. We find an intrinsic spectral slope at 
optical and near-infrared wavelengths of 0.69 ± 0.14 at 0.87 days. The optical 
decay during the first day is almost identical to that of GRB 990123 except that 
GRB 02121 l’s optical afterglow was intrinsically approximately 38 times fainter 
and the transition from the reverse shock to the forward shock may have occurred 
earlier than it did for GRB 990123. We find no evidence for a jet break or the 
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cooling break passing through optical frequencies during the first day after the 
burst. There is weak evidence for a break in the J-band decay between 0.89 
and 1.87 days which may be due to a jet. The optical and infrared data are 
consistent with a relativistic fireball where the shocked electrons are in the slow 
cooling regime and the electron index is 2.3 ± 0.1. The burst appears to have 
occurred in a homogeneous ambient medium. Our analysis suggests that the jet 
of GRB 021211 may have a small opening angle (l?4-4?4) and that the total 
gamma-ray energy is much less than the canonical value of 1.33 x 10 51 erg. If 
this is the case then most of the energy of the burst may be in another form such 
as a frozen magnetic field, in supernova ejecta, or in a second jet component. 
The host galaxy of GRB 021211 is subluminous and has a star formation rate of 
at least 1 A4© yr _1 . 

Subject headings: gamma rays: bursts 


1. Introduction 

The study of gamma-ray bursts (GRBs) took a giant leap forward in 1997 when obser- 
vations with the Dutch-Italian satellite BeppoSAX identified a fading X-ray source inside 
the error box of GRB 970228 (Costa et al. 1997). This rapid localization of a counterpart 
allowed the fading optical afterglow (OA) to be seen (Groot et al. 1997; van Paradijs et al. 
1997). Later that year Metzger et al. (1997) found a redshift of 0.835 for GRB 970508, which 
confirmed the cosmological distances of GRBs. Since then the afterglows of several dozen 
GRBs have been localized. In the cases where the afterglows have been well localized the 
median redshift is approximately two and there is no evidence for extra-Galactic reddening 
of more than approximately one mag along the line of sight to the afterglow. Host galaxies 
have been found for almost every well-localized GRB identified (Bloom et al. 2002). 

The FREGATE, WXM, and SXC instruments on board the High Energy Transient 
Explorer II ( HETE-II ) satellite detected gamma-ray burst (GRB) GRB 021211 (= HETE 
trigger H2493) in the constellation Cancer at 11:18:34.03 UT on 2002 Dec. 11 (Crew et al. 


1 Some of the observations reported here were obtained at the MMT Observatory, a joint facility of the 
Smithsonian Institution and the University of Arizona. Some of the observations reported here were obtained 
with the Magellan 2 (Landon Clay) telescope at the Las Campanas Observatory. Some of the observations 
reported here were taken at the WIYN Observatory, a joint facility of the University of Wisconsin-Madison, 
Indiana University, Yale University, and the National Optical Astronomy Observatories. 
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2002, 2003). The burst consisted of a single peak with a fast rise and exponential decay 
profile. The t $ o durations were 2.30±0.52 seconds at high energies (85-400 keV) and 8.5±0.5 
seconds at low energies (2-10 keV). The energy fluence was (0.96 ± 0.03) x 10' 6 erg cm -2 
between 7 and 30 keV, and (1.98 ± 0.15) x 10~ 6 erg cm -2 between 30 and 400 keV giving a 
hardness ratio of Sx/S y = 0.48 (Crew et al. 2003). The gamma-ray properties of GRB 021211 
were not unusual and fit into the long-soft class of bursts (Kouveliotou et al. 1993). The 
burst was X-ray rich with a spectrum that peaked at Ep^ = 58 keV (Crew et al. 2003). 

The field of GRB 021211 was observed starting 20.74 minutes after the HETE-II trigger 
with the Near Earth Asteroid Tracking Camera on the Oschin 48-inch telescope at the 
Palomar Observatory. These observations led to the announcement of the discovery of the 
OA (Fox et al. 2003). Observations taken even earlier by Wozniak et al. (2002); Li et al. 
(2003), and Park et al. (2002) show that the OA faded rapidly, with a power-law slope of 
« 1.6, between approximately 90 seconds and 11 minutes after the burst. After this it faded 
more slowly, with a power-law slope of 0.82 (Li et al. 2003). Della Valle et al. (2003) found 
a redshift of z = 1.004 ± 0.002 for the host galaxy based on four emission lines from the host 
galaxy. 

Wei (2003) finds that the optical decay of GRB 021211 can be explained if the emission 
before eleven minutes is due to the reverse shock while the emission after that is due to 
the forward shock. They find that the magnetic field strength is unusually low and propose 
that this is why GRB 021211’s OA was faint. However, Kumar & Panaitescu (2003) analyse 
radio and optical data for GRB 021211 and find that the observations before eleven minutes 
after the burst are consistent with a reverse shock if the energy density in the magnetic field 
in the reverse shock is ea r ~ 0.1, higher than usual for a GRB. They also find that the 
local environment around the progenitor is best fit by a homogeneous interstellar medium 
(ISM) with a low particle density (no ~ 10 -2 to 10 -3 cm -3 ). Their model suggests that the 
magnetic field was frozen in during the explosion, and that the magnetic field contains a 
significant amount of the energy from the explosion. 

Della Valle et al. (2003) used deep imaging and spectroscopy to find evidence of a super- 
nova component 27 days after the burst. They find that the optical spectrum is consistent 
with that of SN1998bw and the Type Ic supernova SN1994I. This made GRB 021211 the sec- 
ond GRB — after GRB 030329/SN2003dh (Stanek et al. 2003; Matheson et al. 2003; Hjorth 
et al. 2003) — for which spectroscopic confirmation of a link between GRBs and supernovae 
was obtained. Della Valle et al. (2003) find R = 25.22 ± 0.10 for the host galaxy. 

In this paper we adopt a cosmology with a Hubble parameter of H$ = 70 km s -1 Mpc -1 , 
a matter density of D m = 0.3, and a cosmological constant of = 0.7. For this cosmology 
a redshift of 2 = 1.004 corresponds to a luminosity distance of 6.64 Gpc and a distance 
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modulus of 44.11. One arcsecond corresponds to 16.06 comoving kpc, or 8.02 proper kpc. 
The look back time is 7.75 Gyr. 

2. The Data 

The OA for GRB 021211 is located at R.A. = 08 h 08 m 59^858, Dec. - +06°43'37"52 
(J2000) (Fox et al. 2003), which corresponds to Galactic coordinates of (6 n , Z n ) = (+20?2950, 
215?7486). The reddening maps of Schlegel et al. (1998) give a Galactic reddening of E B ~v ~ 
0.028 ±0.020 mag in this direction. The corresponding Galactic extinctions are Ay = 0.149, 
A b = 0.119, A v = 0.091, A Rc = 0.074, A Ic = 0.053, Aj = 0.026, A H = 0.016, and 
A k = 0.010. 


2.1. Optical Photometry 

Optical imaging data were obtained on 2002 Dec. 12 UT using the Minicam on the 6.5m 
Multi-Mirror Telescope (MMT) at the MMT Observatory. The instrumental gain was 2.7 
e - ADU -1 , the readout noise was 4.5 e - per pixel, and the plate scale was 0"091 per pixel. 

We obtained additional optical imaging data on 2003 Jan. 2 UT using the Low Dispersion 
Survey Spectrograph 2 (LDSS-2) in its imaging mode on the 6.5m Magellan 2 Landon Clay 
telescope at the Las Campanas Observatory. The gain was 3.8 e - ADU -1 , the readout noise 
was 10.2 e“ per pixel, and the plate scale was 07378 per pixel. 

OdL' final set of optical imaging data were obtained on 2003 Jan. 29 UT using Minimosaic 
on the^tVIYN 3.5m telescope at Kitt Peak. The gain was 1.5 e" ADU -1 , the readout noise 
was 5.5 e“ per pixel, and the plate scale was 0 / . / 141 per pixel. 

A log of our observations and the photometry of the OA is presented in Table 1. No 
reddening corrections have been applied to the data in this Table. Figure 1 shows the optical 
data presented in this paper. 

The data were preprocessed using standard techniques for bias and flat-field corrections. 
Photometry was performed by fitting a point-spread function using the DAOPhot II and 
AllStar software (Stetson 1987; Stetson h Harris 1988). The MMT and Magellan data 
were calibrated using secondary standard stars from the catalogue of Henden (2002). Each 
image was calibrated using all of the Henden (2002) secondary standards that could be 
photometered on that image. A weighted mean of the calibrated magnitudes of each of our 
detections of the secondary standards was used to determine the photometric zero point for 



Vega Magnitude 



1 10 100 
Days since 2002 Dec. 11.4712272 


Fig. 1. — These are our optical and infrared data for GRB 021211. The closed circles are 
B-, the closed squares are V-, the closed triangles are R c -, the open circles are J-, the open 
squares are H-, and the open triangles are AVband data. The arrows indicate upper limits. 
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each image. The colour terms in the calibrations were smaller than 5% and did not improve 
the quality of the calibrations. Therefore they were not applied to the photometry. The 
WIYN image from 2003 Jan. 29 was very deep, so the secondary standards of Henden (2002) 
that were in the field of view were saturated. Therefore the zero point for this image was 
determined by matching stars on the i?c-band image obtained with the MMT on 2003 Dec. 
12. 


2.2. Infrared Photometry 

Infrared imaging data were obtained on 2003 Dec. 12, and 13 UT using the Classic 
Camera’s NICMOS3 detector on the 6.5m Magellan 1 (Walter Baade) telescope at the Las 
Campanas Observatory. The gain was 7.5 e~ ADU -1 , the readout noise was 40 e _ per pixel, 
and the plate scale was 0'.'095 per pixel. A dithering sequence with steps of a few arcseconds 
was used for the individual images. The data were preprocessed in the standard way and 
the individual images were coadded to produce a final mean image for each epoch. The 
brightest and faintest ss 5% of the pixels at each location in the stack were rejected in order 
to eliminate cosmic rays and bad pixels that were missed in the preprocessing steps. 

Several faint infrared standard stars from the lists of Persson et al. (1998) were observed 
on each night to determine the photometric zero points of our images. These stars are listed 
in Table 2. The standard star data were preprocessed in exactly the same manner as the 
programme data, and aperture photometry was performed on the individual images. The 
zero-point offset for each standard was computed for each observation of that standard, and 
a weighted mean zero point was computed from all of the standards in each filter on each 
night. We found no evidence for significant time evolution in the zero point over the course 
of a night. Colour terms do not improve the quality of the calibration, so they were not 
used. 

A log of our infrared observations and the photometry of the OA is included in Table 1 . 
No reddening corrections have been applied to this data. Figure 1 shows the infrared data 
presented in this paper. 
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Table 1. Log of the GRB 021211 observations and the results of the photometry. The UT 

date is the middle of each exposure. 


UT Date 

JD - 2450000 

t (days) 

Telescope 

Filter 

Magnitude 

Seeing (") 

Exptime (s) 

Dec. 12.17617 

2620.67617 

0.70494 

Magellan 1/ClassicCam 

K s 

> 21.008 

0744 

660 

Dec. 12.22002 

2620.72002 

0.74879 

Magellan 1/ClassicCam 

J 

21.853 ±0.109 

0748 

3720 

Dec. 12.31163 

2620.81163 

0.84040 

Magellan 1/ClassicCam 

H 

20.921 ± 0.221 

0735 

240 

Dec. 12.31873 

2620.81873 

0.84750 

Magellan 1/ClassicCam 

K s 

> 19.4 

0749 

240 

Dec. 12.33329 

2620.83329 

0.86206 

MMT/Minicam 

V 

23.319 ± 0.150 

0"76 

500 

Dec. 12.34816 

2620.84816 

0.87693 

MMT/Minicam 

B 

23.281 ± 0.225 

0766 

900 

Dec. 12.35336 

2620.85336 

0.88213 

Magellan 1/ClassicCam 

J 

21.700 ± 0.075 

0740 

4680 

Dec. 12.35917 

2620.85917 

0.88794 

MMT /Minicam 

Rc 

22.875 ± 0.086 

0751 

900 

Dec. 13.34512 

2621.84512 

1.87389 

Magellan 1/ClassicCam 

J 

22.812 ± 0.257 

0745 

4320 

Dec. 13.36368 

2621.86368 

1.89245 

Magellan 1/ClassicCam 

J 

> 22.038 

0741 

2400 

Dec. 13.38875 

2621.88875 

1.91752 

Magellan 1/ClassicCam 

K s 

> 19.996 

0742 

210 

Jan. 2.06488 

2641.56488 

21.59365 

Magellan 2/LDSS-2 

Rc 

24.929 ± 0.171 

1713 

300 

Jan. 2.08051 

2641.58051 

21.60928 

Magellan 2/LDSS-2 

V 

25.422 ±0.158 

0775 

300 

Jan. 29.27962 

2668.77962 

48.80839 

WIYN /Mini mosaic 

Rc 

25.070 ± 0.120 

0758 

600 


Table 2. Faint infrared standard stars from the lists of Persson et al. (1998). The 

coordinates are for J2000.0. 


ID 

R.A. 

Dec. 

J 

H 

K s 

SJ9104 

01:03:15.8 

-04:20:44 

11.045 ± 0.005 

10.750 ± 0.005 

10.695 

± 0.010 

SJ9109 

04:18:18.9 

-69:27:35 

11.552 ± 0.002 

11.326 ± 0.002 

11.269 

± 0.002 

SJ9113 

05:10:25.6 

-44:52:46 

11.059 ± 0.005 

10.766 ± 0.005 

10.713 

± 0.005 

SJ9116 

05:42:32.1 

±00:09:04 

11.426 ± 0.009 

11.148 ± 0.009. 

11.058 

± 0.008 

SJ9118 

06:22:43.7 

-00:36:30 

11.723 ± 0.011 

11.357 ± 0.009 

11.261 

± 0.010 

SJ9119 

06:29:29.4 

-59:39:31 

12.114 dfc 0.006 

11.838 ± 0.005 

11.781 

± 0.005 

SJ9121 

06:42:36.5 

-45:09:12 

11.719 ± 0.004 

11.434 ± 0.004 

11.372 

± 0.004 

SJ9134 

08:29:25.1 

±05:56:08 

11.881 ± 0.007 

11.624 ± 0.005 

11.596 

± 0.006 

SJ9137 

09:15:50.5 

-36:32:34 

11.153 ± 0.007 

10.891 ± 0.007 

10.836 

± 0.010 

SJ9187 

23:23:34.4 

-15:21:07 

11.857 ± 0.003 

11.596 ± 0.003 

11.542 

± 0.003 
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3. Results 

3.1. The Optical Light Curve 

To parameterize the optical decay we combined our Ac-band data with the published 
data from Fox et al. (2003), Li et al. (2003), Pandey et al. (2003), and Della Valle et al. 
(2003) (see Figure 2). We fit the data with a double power law of the form 

Rcj&i(t) = -48.77 - 2.5 log 10 (0.5 UU) {{t/t t y ai + (t /t t )~ a2 ] ) + A c , host, (1) 

where t t is the time of the transition from the first to the second power law component, 
f v {t t ) is the flux in the R c band at the time of the transition, R Cj host is the magnitude of 
the host galaxy, and a i: a 2 are the decay slopes before and after the transition respectively. 
This formalism assumes that there are two components to the light from the OA. Both 
components contribute to the total light and each component fades at different rates. We 
fit a double power law instead of a broken power law because the preferred interpretation 
of the early behaviour of the OA of GRB 021211 is a reverse shock. Therefore, the OA 
will contain contributions from both the reverse and forward shocks with the transition time 
indicating when the contribution from each was equal. A broken power law assumes that 
there was a sharp transition between the two decay rates. This is (approximately) valid for 
the transition from a relativistically beamed jet to a non-relativistically beamed jet, but it 
is not valid for the transition between the OA being dominated by a reverse and forward 
shock. 

We fit Eq. 1 to all of the Ac- band data for the first day. We did not include the supernova 
component in our fit since its contribution is negligible at these early times. The magnitude 
of the host galaxy was fixed at Ac, host = 25.22 (Della Valle et al. 2003). Our best fit is shown 
in Figure 2. We find power law slopes of a j = 1.97 ± 0.05 and a 2 = 1.01 ± 0.02; a transition 
time of t t = 5.46 ± 0.80 minutes; and a flux at the transition of f v {t t ) = 0.996 ± 0.214 m.Jy, 
which corresponds to Rc{U) = 16. 12 ±0.23 mag. The fit does not change if we allow the host 
magnitude to be a free parameter. In this case the best fit for Ac, host is 25.25 ± 0.07. Our 
fit is formally not good. The chi-square value is 162.1 for 45 degrees of freedom. Reducing 
this to x 2 /45 = 1 would require that we have underestimated the errors in the magnitudes 
by a factor of approximately two. Since we have used multiple data sets it is possible that 
cross-calibration errors may be present. However, in Sect 3.3 we show that the large residuals 
are present in fits to individual data sets, so the large chi-square value that we find is not 
due cross-calibration errors. 

The double power law does not include the effects of the synchrotron peak frequency, v m , 
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0.001 0.01 0.1 1 10 100 
Days since 2002 Dec. 11.4712272 


Fig. 2. — The upper panel shows our Z?c-band data and the published data for GRB 021211. 
The open circles represent the published data while the solid circles represent our data. The 
line is the best-fitting double power law to the Rc-band data for t < 1 days as described in 
§3.1. The cross at the location of the break shows the size of the la error bars in the fitted 
break time. The lower panel shows the residuals in the fit. The variations during the first 
day are real. 
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moving through the optical. The is expected to happen at early times (t & 5-15 minutes). If 
I'm did pass through the Rc band at approximately the transition time between the forward 
and reverse shocks it will result in the best-fit slopes and transition time in Eq. 1 being 
incorrect. 

Our late-time Magellan data is consistent with there being no variation in the R c - band 
magnitude of the host galaxy between 22 and 49 days after the burst. The only evidence for 
a supernova bump in the optical decay is the photometry of Della Valle et al. (2003), who 
find Rc — 24.48 ± 0.18 at 29 days after the burst, which is approximately 0.5 mag 2.6cr) 
brighter than the combined OA and host galaxy at that time. It should be noted, however, 
that the OA of GRB 030329 did not exhibit a strong rebrightening when SN2003dh was 
near its peak brightness. It is possible that the late-time OA of GRB 021211 behaved in a 
similar manner. This suggests that the lack of an optical bump may not indicator the lack 
of a supernova. 

We find no evidence for a second break in the optical decay during the first day. Our 
J-band photometry suggests that the decay between 21.2 and 45 hours was a = 1.36 4=0.33. 
This is steeper than the decay in the Rc band during the first day and may indicate that a 
break occurred during this time. In § 4.2 we argue that the cooling break is most likely to 
be at X-ray frequencies during this time, so we consider it very unlikely that this apparent 
steepening of the decay is due to the cooling frequency passing through the J band. Therefore 
it is possible that we are seeing the jet break. However, the J band flux increased by ^ 15% 
between 18 and 21.2 hours, so it is also possible that we are seeing a short-term fluctuation 
in the J-band flux (see § 3.3) and not a break in the decay. 


3.2. The Spectral Energy Distribution 

We combined our photometry with data drawn from the literature (Fox et al. 2003; Li 
et al. 2003; Pandey et al. 2003) to construct spectral energy distributions (SEDs) for the 
OA at several epochs. The optical and infrared magnitudes were converted to flux densities 
based on the zero points of Fukugita et al. (1995) and Allen (2000). Each data point was 
corrected for Galactic reddening but not for any reddening that may be present in the host 
galaxy, or in intergalactic space between us and the host. 

The SED was fit by oc x io - 0 - 4 ^), where f v {y) is the flux density at frequency 
v , (3 is the intrinsic spectral index, and A{v) is the extragalactic extinction along the line 
of sight to the burst. The dependence of A(u) on v has been parameterized in terms of 
the extinction in the rest-frame B band, A B . We adopted the three extinction laws given 



-11 - 


by Pei (1992): the Milky Way (MW), the Large Magellanic Cloud (LMC), and the Small 
Magellanic Cloud (SMC). The fit provides (3 and A B simultaneously for each of the assumed 
extinction laws. The unextincted case (A B = 0) was also considered. 

There are enough data to test the four extinction laws only at 0.87 days after the burst. 
We find that our data are consistent with A B = 0 in the host galaxy (x 2 /DOF = 0.5446). 
The three extinction laws of Pei (1992) do not provide a statistically better fit for any value 
of A b in the host. In addition, the lack of curvature in the spectrum after correction for 
Galactic extinction implies that there is very little extinction in the host galaxy along the 
line of sight to GRB 021211. Therefore we believe that there is negligible extinction affecting 
this OA in its host galaxy. 

To test for dust along the line of sight between us and the host we repeated our fits 
allowing the redshift of the dust to be a free parameter. In these cases the best fits for each 
extinction law were for no dust at any redshift in the interval 0 < 2 < 1.004. Therefore, we 
conclude that the dim nature of the OA of GRB 021211 must be intrinsic to the burst itself 
and not the result of absorption along the fine of sight to the burst. 


3.3. Fluctuations in the Decay 

Rapid variability has been seen in the early optical decay of several GRB afterglows. 
GRB 011211 showed fluctuation of approximately 15% over a period of 2.5 horns approx- 
imately twelve hours after the burst (Holland et al. 2002). GRB 021004 exhibited rapid 
variations in flux and colour approximately one day after the burst (Bersier et al. 2003; 
Matheson et al. 2003). Matheson et al. (2003) have also found rapid colour variations in the 
OA of GRB 030329 approximately 0.75 days after the burst. Figure 2 suggests that there 
may be small-scale variability in the R^-band decay of GRB 021211. The anomalous data 
point at 29 days may be due to light from a supernova (Della Valle et al. 2003), so we will 
only consider the optical decay before approximately one day after the burst. The residuals 
about our best-fitting double power law between approximately one minute and one day after 
the burst are not consistent with zero. We find that the y 2 probability of observing these 
residuals about a smooth decay by chance is 2 x 10 -14 . This is sufficient to rule out a smooth 
decay between approximately 1.5 minutes and one day after the burst. A synchrotron break 
passing through the Rc band can produce a bump in the fight curve similar to those seen 
during this time. However, it can not produce the multiple bumps that are observed. We 
can not rule out the possibility that some of the variation seen in the optical decay is due 
to the passage of the synchrotron break. 
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Fig. 3.— The SED of the OA of GRB 021211 on 2002 Dec. 12.341 UT. The filled circles 
represent observed photometry corrected for extinction in the Milky Way. The line represents 
the SED fit assuming no extinction of the host. If we assume that the unextincted spectrum 
follows f„{ u) oc then the best fit has (3 = 0.69 ± 0.14 and Ab = 0 mag. The K s band 
upper limits at approximately 0.87 days are shown as open circles with the time, in days 
since the burst, given in parentheses. Note that the predicted A s -band flux at 0.87 days 
(5.367 ^Jy) is greater than the upper limit on the AT s -band flux at 0.70 days (2.638 ^Jy) 
which implies that the infrared afterglow became brighter between 0.7 and 0.87 days. 
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To test that this result is not due to a single data set, which would suggest that the 
fluctuations are an artifact of the data reduction, we computed the \ 2 probability of observing 
the residuals about a smooth decay by chance for each data set. The Fox et al. (2003) data 
have a probability that the residuals were due to chance of 3 x 10~ 9 . The Li et al. (2003) 
data have a probability that the residuals were due to chance of 1 x 10 -8 . The Pandey et al. 
(2003) data are consistent with random fluctuations about a smooth decay. It is not possible 
to determine if the systematic offset of +0.14 mag from the best-fitting smooth decay seen 
in the Pandey et al. (2003) data is real or an artifact of their data reduction. The offset is 
consistent with the fluctuations seen in the other data sets, so we have no reason to believe 
that it is not real. 

The upper limit on the AT,- band flux at 0.70 days is 2.638 /iJy. This is 2.729 p.Jy (0.77 
mag) fainter than the flux predicted by extending the fitted spectrum at 0.87 days to the 
K s band. This suggests that the infrared afterglow has brightened during this time. The 
J-band flux also appears to have brightened during this time. If the increase in the J-band 
flux is part of an achromatic fluctuation in the flux then it predicts K s = 20.43 at 0.705 
days. However, the upper limit is K s > 20.998 at this time. This suggests that the infrared 
fluctuations are not achromatic. 

Pandey et al. (2003) report spectral slopes, corrected for Galactic reddening, at three 
epochs between 0.13 and 0.86 days after the burst. Their slopes are consistent with a constant 
spectral slope dining this time of (3 = 0.56 + 0.13 (B—V — 0.13) at the 93% confidence level. 
At t = 0.1 days Fox et al. (2003) find /? = 0.98 (B — V = 0.23) after correcting for Galactic 
reddening. We find an intrinsic spectral slope of 0 — 0.69 ±0.14 (B — V = 0.28) at t = 0.87 
days. The only evidence for a change in colour are the data of Fox et al. (2003) who derived 
their spectral slope from B- and A', -band photometry. When we combine our colour data 
with that in the literature (Fox et al. 2003; Li et al. 2003; Pandey et al. 2003) we find no 
evidence for a change in colour between 0.1 and 1.0 days. However, we note that there is 
too little published data to be able to make any strong statements about colour variations 
during this time. 


4. Physics of the Burst 
4.1. Energy Considerations 

We computed the isotropic equivalent energy of GRB 021211 using the observed HETE- 
II fluences listed in Crew et al. (2003). We applied a cosmological K correction (Bloom 
et al. 2001) to the energy in each band in order to correct it to the 20-2000 keV band then 
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averaged the six results. The resulting K - corrected isotropic equivalent energy is E lso = 
(1.0 db 0.1) x 10 52 erg. 

We find no evidence for a break in the optical decay due to beaming before approximately 
one day after the burst. However, Kumar & Panaitescu (2003) note that the radio flux density 
at 8.5 GHz radio flux at 10 days is a factor of 5.4 lower than expected. This discrepancy 
can be explained if a jet break occurred between one and ten days after the burst. If the 
steepening of the J-band decay that we observe between 21.2 and 45 hours after the burst 
is due to a jet break then the time of the break is constrained to 0.89 and 1.87 days. 

The opening angle of a GRB jet can be estimated using Rhoads (1999) and Sari et al. 
(1999). This requires a knowledge of the circumburst density, which is not well understood 
for GRB 021211. Detailed modeling by Kumar h Panaitescu (2003) implies 10~ 3 < n < 10 -2 
g cm -3 . If we assume an efficiency of r/ 7 = 0.2 and a jet break between one and ten days 
the opening half-angle of the jet is 1?4 < Qj < 4?4. If the Kumar &; Panaitescu (2003) 
density range is correct then GRB 02121 l’s jet has one of the smallest known half-opening 
angles. The beaming-corrected energies that corresponds to the estimated range of jet angles 
is 3 x 10 48 erg < E y < 3 x 10 49 erg, which is a factor of approximately 40 to 400 smaller than 
the canonical value of Bloom et al. (2003). In order for GRB 021211 to have the canonical 
energy the jet half-opening angle needs to be & 30°. If this is correct then either the jet break 
will not occur until t ^ 10 days after the burst or the ambient density exceeds ^ 48 000 
cm -3 . Both of these scenarios are unlikely so we believe that GRB 021211 may have been 
underluminous at gamma-ray energies. 


4.2* The Cooling Frequency and Local Environment of GRB 021211 

Most GRBs are well fit by relativistic fireball models with an electron index of p ~ 
2. 3-2. 5 (van Paradijs et al. 2000). Panaitescu & Kumar (2002) present models for ten 
GRBs and find a mean electron index of p = 1.9 with five GRBs having p < 2. This is 
problematic since electron indices of less than two represent infinite energy in the standard 
relativistic fireball model (Meszaros 2002). This problem can be avoided by introducing an 
upper limit for the electron energy distribution (Dai h Cheng 2001). Detailed modeling 
of the acceleration of particles in highly relativistic shocks predict an electron index of 
approximately 2.3 (Achterberg et al. 2001). 

We estimate the electron index for GRB 021211 using the slope of the optical decay, 
and the spectral index, of GRB 021211 at 0.87 days. The optical decay has a slope of 
a = 1.01 ±0.02 and the dereddened slope of the SED at optical wavelengths is (3 = 0.69±0.14 
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(B — V = 0.28 ± 0.03). We used the relationships of Sari et al. (1998) (for a homogeneous 
interstellar medium (ISM)) and Chevalier &; Li (1999) (for a pre-existing stellar wind) to 
predict p at this time. These predictions are listed in Table 3 for several stages of the 
evolution of a synchrotron emission spectrum. 

Cases 5, 6, 11, and 12 can be ruled out because they predict a spectrum with an opposite 
slope to what is observed. Cases 9 and 10 can be ruled out because they predict that the 
optical decay has a slope of a = 0.25, which is not consistent with the observed slope. The 
SED at t = 0.87 days is barely consistent (« 2.5 a) with a prediction of 0 = 1.01 (cases 1, 
2, 7, and 8) and 0 = 0.34 (case 4). It is, however in good agreement with 0 = 0.67 (case 
3). This suggests that GRB 021211 was in the slow cooling regime with the synchrotron 
frequency, i/ m , below the optical and the cooling frequency, v c , above the optical 0.87 days 
after the burst. 

This conclusion is supported by the detailed modelling of Kumar k. Panaitescu (2003). 
They found that the early- and late-time optical data is consistent with a low-density, ho- 
mogeneous external medium and inconsistent with a wind-stratified external medium. Our 
broad band spectral modelling supports this conclusion. However, we are unable to rule out 
the possibility that the burst expanded into a shocked wind, or into an environment where 
the equatorial wind is significantly stronger than the polar wind. 

If the burst did expand into a homogeneous medium then the observed rate of decay of 
the OA at t = 0.87 days (a = 1.01 ± 0.02) and case 3 of Table 3 implies that the electron 
index is 2.3±0.1 and the intrinsic spectral slope at this time is 0.67±0.01. This is consistent 
with the dereddened spectral slope at this time (0 = 0.69 ± 0.14) that we determined in 
§ 3.2. Eq. 11 of Sari et al. (1998) implies that the cooling frequency at t = 0.87 days in a 
homogeneous medium is v c — 2.9 x 10 12 €g 3/,2 n _1 Hz. For the case of a homogeneous ISM 
Kumar k Panaitescu (2003) find es f « 0.001 in the forward shock at 11 minutes and an 
environmental density of 10 -3 < n < 10 -2 g cm -3 . Using these value we predict that the 
cooling frequency at t = 0.87 days should be approximately 9 x 10 18 ^ c <9x 10 19 Hz, 
which is in the high-energy regime (approximately 37-370 keV). However, the modeling of 
Kumar k Panaitescu (2003) suggests that €b } decreases by a factor of two between eleven 
minutes and ten days after the burst, so locating the location of the cooling frequency is 
somewhat uncertain. Regardless, the Kumar k Panaitescu (2003) model is inconsistent with 
v c being near the optical at this time. Therefore, we do not believe that the cooling break 
moved through optical frequencies between 0.1 and 1.0 days after the burst. 

The emerging picture of the evolution of this burst is a transition from the emission 
being dominated by the reverse shock to emission being dominated by the forward shock 
at approximately five minutes after the burst. The reverse shock has an optical decay of 
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ai = 1.97 ±0.05, the forward shock has an optical decay of = 1.01 ±0.02 and the spectral 
index of the forward shock photons is j3 = 0.69 ±0.14. The spectral and decay slopes 
imply that the electron index is p = 2.3 ±0.1 and are consistent with the local environment 
dominated by a homogeneous ISM. 


5. Discussion 

5.1. A Comparison with GRB 990123 

The overall shape of the optical decay of GRB 021211 appears to be very similar to that 
of GRB 990123 (e.g., Holland et al. 2000). To quantify this we fit Eq. 1 to the early-time Rc- 
band data from Fruchter et al. (1999) and the data used by Holland et al. (2000). The best 
fit occurs with power law slopes of Qi = 2.25 ± 0.48 and a 2 = 1.15 ± 0.05, a transition time 
of t t = 1.3 ± 1.0 minutes, and an f?c-band magnitude at transition of Rc{t t ) — 10.91 ± 1.33. 
The goodness of fit was X 22 = 0.4335. The magnitude of the host was fixed since we are 
interested in determining the slopes and transition time under the assumption that there are 
two mechanisms (the forward and reverse shocks) contributing to the flux, not the magnitude 
of the host. 

The OA of GRB 990123 made the transition from being dominated by the reverse shock 
to being dominated by the forward shock faster than the OA of GRB 021211 did. The rest- 
frame transition times are 2.7 minutes for GRB 021211 and 0.5 minutes for GRB 990123. 
However, this difference may be nothing more than an artifact of the data. The GRB 990123 
data consists of 26 data points from several sources, so there may be a systematic offsets in 
the photometric calibrations. If we force the GRB 990123 fit to have a rest-frame transition 
time of 2.7 minutes (the same as GRB 021211 does) then the best fit double power law has 
ai = 1.79 ± 0.10 and <22 = 1.06 ± 0.02 with a goodness of fit of X 22 = 0.5143. This is only 
marginally worse than the best fit with the transition time is a free parameter and the decay 
slopes are only 1.7a different from the fit to the GRB 021211 data. This suggests that the 
overall temporal evolution of the optical decays of both bursts was the same to within the 
limitations of the data. 

These two OAs differ in that GRB 021211 shows strong evidence for rapid variations 
about a power law decay during the first day while GRB 990123 exhibits no variations 
(Holland et al. 2000). This may be an artifact of the sampling rate of the GRB 990123 
data, or cross-calibration issues. However, f the smooth optical decay of GRB 990123 is real 
then it may indicate that the circumburst medium around GRB 990123 is smoother than 
that around GRB 021211 (Wang h Loeb 2000). Alternately, GRB 021211 may have had 
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additional energy added via refreshed shocks (Panaitescu et al. 1998) while GRB 990123 did 
not. 

Figure 4 shows the Rc-band data for GRB 990123 with the best-fitting double power 
law superimposed. The poor temporal coverage of the data makes it difficult to determine 
the transition time. This Figure also shows the Rc-band data for GRB 021211 with the 
fit for GRB 990123 superimposed. With the exception of the short-timescale fluctuations 
the shape of the GRB 990123 fit is in good agreement the optical decay of GRB 021211. 
The data for GRB 990123 has been adjusted to correct for the difference in redshift between 
the two GRBs. The fluxes with shifted to match the observed Rc band of GRB 021211 by 
assuming a power law spectrum with (3 = 0.75 (Holland et al. 2000) for GRB 990123. Next 
the difference in distance modulus (Ap = 1.25 mag) was applied. Finally, a K correction 
was applied to compensate for the difference in redshift between the two bursts. After these 
effects had been taken into account the OA of GRB 021211 was 3.95 mag fainter than the 
OA of GRB 990123, which corresponds to a difference in intrinsic luminosity of 38 times. 
Except for the small-timescale fluctuations this difference is constant between approximately 
90 seconds and one day after the burst. 

The spectra of these two OAs are similar. Holland et al. (2000) used broadband pho- 
tometry to find / 3 = 0.75 ± 0.07 for GRB 990123 between zero and three days after the 
burst. This is almost identical to the value for GRB 021211. Andersen et al. (1999) found 
(3 = 0.69 ± 0.10 from a spectrum taken 48 minutes after the burst. The similarity in the 
spectra of the two bursts implies that they had similar electron indices. 


5.2. The Host Galaxy 

The location of the OA in our 2003 Jan. 29 UT W1YN image was determined by a 
comparison with HST /ACS/WFC images taken on 2003 Dec. 24/25 UT (Fruchter et al. 
2002). Our late- time image has a point source with Rc = 25.07 ± 0.12 at the location of the 
OA. Fruchter et al. (2002) find that this source is resolved with an intrinsic full width at half 
maximum of ~ O'.' 13 (~ 1 proper kpc) once the point-spread function of the ACS is removed. 
This is beyond the ability of WIYN to resolve so it appears as a point source in our data. 
The HST/ ACS/WFC data give an AB magnitude in the F606W filter of 25.3 ± 0.2, which 
corresponds to Rc ~ 25.1. We find Rc = 24.929 ± 0.171 at t — 21.6 days after the burst 
and Rc = 25.070 ± 0.120 at 48.8 days. 

The V — R colour of the afterglow at t = 21.6 days is 0.47 ± 0.23, after correcting for 
Galactic reddening, which corresponds to a spectral slope of (3 = — 1.4 ± 1.2. The data in 



AB Magnitude 



Fig. 4. — This Figure shows the GRB 990123 R c -band data (open squares) with the best- 
fitting double power law as described in § 5.1. The published GRB 021211 Re- band data 
(open circles), and our Re- band data (closed circles), is also shown. The best- fitting double 
power law to the GRB 990123 Rc-band data superimposed has been shifted 3.95 mag fainter 
so that it overlies the GRB 021211 data. The AB magnitude system is used in this Figure. 
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de Vaucouleurs et al. (1991) suggest that the spectral slope of the Small Magellanic Cloud 
(SMC) at near- ultraviolet wavelengths is /?smc ~ —1.6. This is consistent with the observed 
colour of the OA at similar rest-frame wavelengths 21.6 days after the burst. 

However, Della Valle et al. (2003) showed that the OA has a supernova component 
that reached its maximum light at ~ 25 days after the burst. They find that, at this 
time, the host and the supernova contribute approximately equally to the light while the 
OA contributes less than approximately 5% of the light. Their best estimate of the host 
magnitude is Rc = 25.22 ± 0.10. The corresponding absolute rest-frame magnitude of the 
host is My ~ —19. 

Lilly et al. (1995) find a typical magnitude of (M B ) AB = —21.4 for blue galaxies at 
0.75 < z < 1.00 if ( //(). f2 m , Sl\) = (50, 1,0). For our adopted cosmology this corresponds 
to (M B ) AB = —22.2. If we assume that the host galaxy has a power-law spectrum with 
/3 g ai = —1.4 (as implied by the V—R colour of the afterglow at t = 21.6 days) then the rest- 
frame 5-band luminosity is L B ~ 0.1 L* B where L* B is the rest-frame 5-band luminosity of a 
typical blue galaxy at 2 m 1. Therefore the host of GRB 021211 appears to be subluminous. 
This is consistent with Hogg k Fruchter (1999) who found that GRB host galaxies tend 
to be subluminous. Our comparison of the total 5-band luminosity of the host to L* B is 
somewhat uncertain since M* is highly correlated with the slope of the faint end of the 
galaxy luminosity function as well as with its normalization (Lilly et al. 1995). 

At a redshift of z = 1.004 a rest-frame wavelength of 2800 A corresponds to the ob- 
server’s V band. Therefore we are able to use Eq. 2 of (Madau et al. 1998) to estimate the 
integrated star formation rate in host galaxy from the rest-frame continuum flux at 2800 A. 
We find a star- formation rate of « 1 _M 0 yr -1 . This calculation depends on the details of the 
initial mass function, and it assumes that there is no extinction in the host galaxy. We know 
from the observed SED of the OA (see Sect. 3.2) that there is negligible extinction in the 
host along the line of sight to the burst. However we do not know what level of extinction 
is present in the rest of the host. Therefore we estimate that the integrated star-formation 
rate in the host galaxy of GRB 021211 is > 1 MqYX^ 1 . 

The specific star-formation rate per unit blue luminosity of the host galaxy of GRB 021211, 
if we ignore extinction within the host, is & 10 Al 0 yr 1 , which is similar to that of 
other GRB host galaxies. Our imaging data do not have high enough resolution to deter- 
mine morphological properties of the host. That will require HST imaging of the host after 
the OA has faded. However the early HST/ ACS imaging of Fruchter et al. (2002) suggests 
that the host has a diameter of « 1 proper kpc. This implies a star-formation rate density 
of > 1.3A4 0 yr” 1 kpc -2 in the host. This is at the low end of the range of values (1-1000 
A40yr -1 kpc -1 ) that characterize a starburst galaxy (Kennicutt 1998). 
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6. Conclusions 

We present BVRcJHK s photometry of the OA of GRB021211 taken at the Magellan, 
MMT, and WIYN observatories. These data were taken between 0.7 and 50 days after the 
burst. The broad-band optical/infrared SED yields an intrinsic spectral slope of 0.69 ± 0.14 
at 0.87 days after the burst and we find no evidence for colour evolution in the optical 
between 0.1 and 1.0 days. There is weak evidence for brightness fluctuations in the J and 
K s bands approximately 0.7 days after the burst. We find no evidence for extragalactic 
extinction along the line of sight to the burst. The optical light is dominated by power law 
decay with an index of 1.97 until 5.46 minutes after the burst. After this time the decay is 
dominated by a power law with an index of 1.01. This is consistent with the early time flux 
being dominated by a reverse shock and a transition to domination by the forward shock 
5.46 minutes after the burst. 

There is evidence for rapid fluctuations in the flux about the smooth power-law decay 
similar to what has been seen in GRB 011211 (Holland et al. 2002), GRB 021004 (Bersier 
et al. 2003; Matheson et al. 2003), and GRB 030329 (Matheson et al. 2003). These fluctua- 
tions may also include colour variations, but there is insufficient data to determine this is a 
convincing way. 

The spectral slope has been combined with the observed Ac-band optical decay to 
determine that the shocked electrons are probably in the slow cooling regime with an electron 
index of 2.3 ±0.1, and that the burst probably occurred in a homogeneous ISM. There is 
weak evidence for a jet break between 0.89 and 1.87 days after the burst, and (Kumar & 
Panaitescu 2003) suggest that the jet break may have occurred before ten days based on the 
radio flux at 8.5 GHz. If the jet break did occur before ten days then the half-opening angle of 
the jet is 1?4 < 9j < 4?4, which is one of the smallest opening angles of any GRB jet (Bloom 
et al. 2003). These angles imply that the total gamma-ray energy in the burst was 3 x 10 48 
erg < E~f < 3 x 10 49 erg, which is much less than the canonical value of Bloom et al. (2003) 
yet similar to that inferred from the first jet break in GRB 030329. This suggests that most 
of the energy in GRB 021211 may be not be in the OA, but in a different component such 
as a frozen in magnetic field (cf Kumar Sz Panaitescu 2003) , kinetic energy in the supernova 
ejecta, or a second jet component. 

The optical decay of GRB 021211 during the first day after the burst was very similar 
to that of GRB 90123, only the OA of GRB 021211 was instrinsically ~ 38 times fainter. 
Both OAs had an initial steep decline with a power law index of ~ 2. This is consistent with 
the OA being dominated by emission from a reverse during the first few minutes after the 
burst. The OA of GRB 990123 made the transition from being dominated by the light from 
the reverse shock to being dominated by the light from the forward shock 1.3 ± 1.0 minutes 
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after the burst (0.5 minutes in the rest frame of GRB 990123). The transition occurred at 
5.46 ± 0.80 minutes for GRB 021211 (2.7 minutes in the rest frame). After the transition 
each OA is dominated by forward shock emission with decay indices of m 1. The similarity 
in the decays of these two OAs suggests that the environments and physics of GRB 990123 
and GRB 021211 were similar. 

The rapid localization of GRB 021211 and the near-continuous monitoring of its OA 
during the first day after the burst occurred has allowed the early time evolution of this burst 
to be studied in detail. This detailed coverage of the optical decay has made it possible to 
recognize the similarity between GRB 021211 and GRB 990123. The OA of this burst faded 
very rapidly for the first several minutes, and had the initial afterglow searches been delayed 
by even one hour the OA would have been « 5 mag fainter. It is questionable whether or not 
the OA would have been discovered in that case. If the OA had not been detected rapidly 
this would have been classified as a dark burst. Further, without the prompt discovery of 
the OA the rapid early time decay would not have been identified. Unfortunately there is 
very little data for the period between one and 20 days after the burst, so it is not possible to 
determine when the jet break occurred. Daily monitoring of OAs during the first few weeks 
will be possible using the the Ultraviolet/Optical Telescope on the upcoming Swift mis sion. 
This will monitoring will make it much less likely that critical periods of the decay, such as 
the jet break, will be missed for bursts that occur during the Swift era. 
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Table 3. Predicted electron and spectral indices between 12.1 minutes and 21.3 hours 
after the burst assuming a = 1.01 ± 0.02. The electron index is predicted from the 
observed decay. The spectral index is derived from the predicted value of p. 


Case 

Model 

Environment p 

0 

Comments 

1 

Vm < v c < V 

ISM 

2.0 ±0.1 

1.01 ±0.01 

2 


Wind 

2.0 ±0.1 

1.01 ±0.01 

3 

l 'm<V <V C 

ISM 

2.3 ±0.1 

0.67 ±0.01 

4 


Wind 

1.7 ±0.1 

0.34 ± 0.01 

5 

V < Vm < Vc 

ISM 


-0.33 

Rising spectrum 

6 


Wind 


-0.33 

Rising spectrum 

7 

Vc < Vm <V 

ISM 

2.0 ±0.1 

1.01 ±0.01 

8 


Wind 

2.0 •0.1 

1.01 ±0.01 

9 

V C < V < Vm 

ISM 


0.5 

a = 0.25 

10 


Wind 


0.5 

a = 0.25 

11 

V <V C <Vyn 

ISM 


-0.33 

Rising spectrum 

12 


Wind 


-0.33 

Rising spectrum 


